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We present a theory to model carrier mediated ferromag- 
netism in concentrated or diluted local moment systems. 
The electronic subsystem of the Kondo lattice model is 
described by a combined equation of motion / coherent 
potential approximation method. Doing this we can cal- 
culate the free energy of the system and its minimum ac- 
cording to the magnetization of the local moments. Thus 
also the Curie temperature can be determined and its de- 
pendence on important model parameters. We get qual- 
itative agreement with the Curie temperatures' experi- 
mental values of Gai_2;Mna;As for a proper hole density. 



1 Introduction Ferromagnetism is a solid state phe- 
nomenon which arises from the interplay of magnetically 
active atoms or their electrons, respectively. Naively one 
would think that high Curie temperatures are connected to 
a high density of those atoms. Thus it was suprising when 
relatively high Curie temperatures |Il]|2]|3] in diluted mag- 
netic semiconductors (DMS) with Tcs over lOOK were 
discovered - even at small amounts of magnetic atoms. 
This finding and the high spin polarization of the carriers 
made the DMS a possible candidate for spintronic devices ID . 
It is a major task to understand these materials and thereby 
improve them to eventually achieve room temperature mag- 
netism. 

One prominent DMS is Gai_j;Mn^As which has a Tc up 
to w 173K for a concentration of x = 9% of manganese 
ionsQ. Manganese, which replaces dominantly gallium, 
has a partly filled 3d^-shell, in contrast to gallium's com- 
pletely filled 3(i^°-shell. To calculate the effect of this dop- 
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ing several methods can be used. Ab-initio calculations are 
able to calculate the full density states of a specific ma- 
terial. This gives insight into the position of the impurity 
levels within the bulk material. Additionally Curie tem- 
peratures of different DMS can be calculated accurately ||6] 
[Tj. Besides ab-initio methods, model studies are useful to 
find the main mechanism which leads to ferromagnetism 
in DMS. For Gai_2;Mna;As the Sd-states of manganese lie 
deep below the valence band||4l|6l of the host material and 
are therefore assumed to lead to a localized spins Si at the 
impurity site R; with spin quantum number S = 5/2. Al- 
though the states are deep below the Fermi level they can 
interact with carriers (holes) in the valence band which in- 
fluences its band structure. This results in an effective ex- 
change between the localized spins Si (even though they 
are far away from each other at strong dilutions) which can 
lead to ferromagnetic order 

One of the major problems is to have control over the car- 
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rier density (n for electrons or p for holes). Nominally one 
manganese ion provides one hole but due to compensation 
effects the ratio 7 ~ p/x is usually smaller than one. As 
the mean-field Zener model predicts Tc ^ p^ this seems 
to reduce the Curie temperature. In contradiction experi- 
ments show|[8][9l that a hole compensation is rather neces- 
sary for a finite Tq- The same is true for theoretical ap- 
proaches beyond the mean-field approximation O 1 Oil 1 1 |[T2l . 
To calculate Tc the Kondo lattice model (KLM) is often 
mapped onto a Heisenberg model, which is then solved 
with different approximation methods. Contrary to that we 
propose a procedure how to calculate the free energy ex- 
clusively from the electronic subsystem so that we get the 
Curie temperature directly from the KLM. 
The paper is organized as follows. First we introduce the 
model Hamiltonian. After that we describe how we calcu- 
late the internal energy of the system by an equation of mo- 
tion approach. In the case of a diluted system we combine 
that with a coherent potential approximation. As the next 
step we use the internal energy to get the free energy where 
we also need the entropy at T = OK. From the free energy 
we calculate Curie temperatures in dependence of several 
model parameters. These are compared with experimental 
data. 

2 Model We use the Kondo lattice model to describe 
the interactions between local moments Si and itinerant 
conduction electrons with spin cTi, with the standard Hamil- 
tonian 

n =Ho + Hsf 

= ^»J'="tCj<t - S, • CT; . (I) 

The itinerant electrons, represented by the creators (anni- 
hilators) c^^K move from lattice site to Rj due to the 
hopping Tij . They prefer to align parallel (J > 0) or an- 
tiparallel (J < 0) to the local moment on the same site. 
In the context of DMS the localized moments are formed 
by the d-states of the magnetic impurities while the itiner- 
ant carriers move in the valence band (holes) or in the con- 
duction band (electrons). In Gai_j;Mna;As the hole created 
by the replacement of Ga with Mn is within the three an- 
tibonding sp — d-levels with dominant As 4p character||4|. 
This actually requires a multiband treatment but to simplify 
the calculations the valence band is approximated in a one- 
band model. The conduction band is completely neglected. 
Since the magnetic impurities are randomly distributed within 
the host material we have to include disorder to our model. 
This is done via an introduction of a random variable Xi = 
0, 1 indicating if a magnetic impurity is present at lattice 
site Ri. The concentration x of the magnetic atoms is fixed 
by the requirement x ^ Xi. Besides the magnetism 

of the impurities there can be other effects on the electrons lfT3l 
[Tn . e.g. different atomic levels for magnetic (a =M) 
and non-magnetic sites (a =NM). We choose = T}f = 



while Tq = Tfj^ may be unequal zero. We end up with 
the Hamiltonian 

- ^ E {^'y^t^^'y + S^4aC^s) (2) 

where we wrote the scalar product in dO with the help of 
Pauh matrices a and <Ti = ^ 12cr(j' cf^^aa'Cia' explicitly. 
There is now an Ising term z^S^rii^ (za = ±1, riia- ~ 
cfaCia) and a spinflip part 5f c+ (Sf = Sf + z^iSf, 
a = —a). 

Since our model is not restricted to small doping x we are 
able to calculate the whole range a; = ... 1 to show gen- 
eral trends of the influence of dilution. 

3 Electronic part In this work the key quantity is the 
internal energy U = (H). It can be derived in the KLM 
from the single-electron Green's function (SE-GF) 

= -z©(<-<')([c..W,c+(i')]±) (3) 

with the thermodynamic average (...) and the usual step 
function 0{t — t'). Operators are in the time dependent 
Heisenberg representation. The GF can be Fourier-transformed 
to its energy representation 

G,,,{E) =jd{t- t')G,,„{t - i')e'^^*-*'^ (4) 

=((cm;c+,)> • 

In the KLM the internal energy then is simply given as 

^ = - ]^ E J dEEf_ {E)lmGua {E) , (5) 

-oo 

where the Fermi function f-{E) ~ (^el^i^-t^) -\- ij^i con- 
tains the chemical potential ^ and the inverse temperature 
(3^^ — ksT. Thus, to calculate U, we need to find an ex- 
pression for Gija (E). In a first step this will be done for the 
concentrated lattice {xi = 1, Vi) and then will be extended 
to diluted lattices. 

3.1 Single electron Green's function in the con- 
centrated lattice To get the SE-GF we will perform a 
moment conserving decoupling approach (MCDA lfT4l ) of 
higher Green's functions appearing in several equations of 
motion (EOM). Let us start with the EOM of G^j^, which 
reads 

E (ES^i - Tu) Gi.AE) = % + (([c,,, c+ )) . 

I 

(6) 
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We can introduce a self-energy 
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^+ W — 



= Y,M,UE)Gij.{E) (7) 
I 



which leads, after Fourier transformation, formally to the 
solution 



G'k<x(-B) = 



1 



^-e(k)-Afka(S) 



(8) 



The commutator in ^ can also be calculated expUcitly, 
yielding 

J 



-(z,(5^)G,j,(i?) + (9) 
r.kjAE) =(( {Sf - (St)) c,,; c+ )) (10) 



(11) 



with the reduced Ising function Fikja (E) and the spinflip 
(SF) function Fikja{E). These functions can also be cal- 
culated by their EOMs. Within these equations we treat di- 
agonal (i = k) and non-diagonal terms (i ^ k) differently. 
For the non-diagonal elements we replace approximately 

h±a, H,f]_ ^ Ma±a{E)ci±„ (12) 



as in Eq. (|7]i, which is only strict for Gij±a{E), of course. 
As an example: 



I 



Mkia{{S!ciy,c+J . (13) 



The higher Green's function is therewith traced back to a 
linear combination of Ising functions. It is easy to see that 
the functions on the left- and the right-hand side have the 
same pole structure differing only by the respective spec- 
tral weights. The latter is approximately compensated by a 
self-consistent determination of the expansion coefficients 
which are just the elements of the central self-energy ma- 
trix. Higher Green's functions with commutators contain- 
ing only spin operators (e.g. [Sf , Hsf]_), are neglected, 
since they scale with magnon energies which are much 
lower than those of the electrons [ 14 1. The commutators 
of the diagonal terms, however, are calculated explicitly 
which leads to further Green's functions {{AiaCi±a', Cjg-))' 
where Ai^ stands for different combinations of spin and 
electron operators. These functions can be expressed in 
some limiting cases (e.g. 5 = |, {S') = 5, n = or 
n = 2) rigorously by the three Green's functions men- 
tioned above suggesting the general ansatz 

((A,crC.,±cr;C+^)) ^aAaGijcriE) + 

+ bAar,^ja{E) + CAaFujaiE) (14) 



which we approximate to be valid for all parameter constel- 
lations. The prefactors UAa, ■ ■ ■ are fixed by the require- 
ment of conservation of the first spectral moments of the 
Green's functions and their according expectation values. 
Since the moments are conserved we then get a correct 
high energy behavior of these functions ill 511 161 . 
After performing the mentioned approximations we get unique 
solutions for the Ising- and spinflip Green's functions (de- 
tails in Ref. lfT4ll ). These can be written in the form 

XujAE) ^ CxAE)G,,^iE) ,X = r,F (15) 

and we see that they are proportional to the one-electron 
Green's function. This allows us, by comparing Eqs. (I6]|2J 
and (|9]l, to identify the prefactors Cxa{E) in (flST l with the 
self-energy 

M^{E)^~'-^z,{S')-^{z„Cra{E) + CFa{E)) , 



j^jmcLin- field 



(16) 



which splits into a mean-field and a more complicated "many- 
body" part. The latter is defined by the Cxa{E), which 
depend on several expection values. These either contain 
only spin operators or also electronic operators. The "pure 
spin operator" values (e.g. ((5^)^)) can all be expressed 
by the well known formula of Callen lTTl as functions of 
the local moment magnetization {S^). To do this we im- 
plicitly assume that a mapping of the KLM onto a Heisen- 
berg model is in principle possible, but we do not have 
to perform it explicitly. All other expectation values (e.g. 

(5f c^^Cicr), {S^Ui^)) can be calculated with the 
spectral theorem from the known Green's functions Gija-, 
Fiija and Fuj^. But the Cxa{E) do not only depend on 
scalars but also on the local Green's function G±cr{E) = 
1 /N Gk±a{E) and the self-energy M±a{E) itself. These 
dependencies necessitate a self-consistent calculation of the 
self-energy. 

Especially the dependence on the local Green's function 



M,{E)=AUE,G±^{E)) 



(17) 



will play a crucial role when we consider diluted systems. 
Since we want to focus on the interaction of local moments 
and carriers, we have neglected a direct influence of the 
Coulomb repulsion in the conduction band. An addition 
of that would be especially necessary to describe narrow 
band materials like the manganites[ 18|. The inclusion of 
a Hubbard-type Coulomb term could lead to extra effects 
in magnetic ordering (band-magnetism), which are not ac- 
counted for in this work. 

3.2 Single electron Green's function in the diluted 
lattice The diluted KLM can be described as a binary al- 
loy. The components are defined by the sites with or with- 
out magnetic impurities. In the coherent potential approx- 
imation (CPA IIT9I ) the configurationally averaged Green's 
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function of the total system is given by 



equations is closed. The self-consistency cycleQ now reads 



(18) 



The CPA self-energy Scr (E) can be derived self-consistently 
from 







'la 



l-('7; 



NM 



(19) 



As the concentration of the magnetic impurities x is a pa- 
rameter we only need to determine the potentials t^^^nm 



to solve the problem. We set rj^ 



NM 



T]^^ for the non- 



magnetic sites. For the magnetic sites we choose the self- 
energy (fT6] l as an energy dependent potential rf^{E) = 
M„{E) (dynamical alloy analogy, DAA 12011 101 ). In Sec- 
tion l3.1l we calculated the self-energy Ma{E), put it into 
formula (O to get G\^a{E) and calculate Ma{E) again un- 
til self-consistency is achieved. But we have to keep in 
mind that the self-energy was derived for a concentrated 
lattice. In a diluted system the magnetic component is em- 
bedded in the total system. Thus we change the self-consistency 
cycle a bit. 

Since the main change of the Hamiltonian (|2| for the di- 
luted system is only the inclusion of scalar random vari- 
ables Xi, we assume that the equations of motion of the 
magnetic subsystem have the same formal structure as in 
Section im That means especially that the self-energy M^r {E) 
is a function of the local Green's function G±a{E) (Eq. 
(fTTli). As we want to calculate only the self-energy of the 
magnetic subsystem it should only depend on the local 
Green's function of this subsystem. Thus we replace G±a [E] — 
{E) in Eq. (HTJl yielding 



S,{E) ^ {G^a{E)) 
^ G'^\E) 



Ma{E)=v^ 

Sa{E)... , 



where the self-energy determination and the CPA formal- 
ism are strongly correlated. 

4 Free energy In the previous sections we determined 
the Green's function and the internal energy, but only for a 
given magnetization (5*"). To decide which magnetization 
the system actually favors we want to calculate the free 
energy and find its minimum, according to (S*^). To ab- 
breviate the notation and avoid confusion with the entropy 
Sm{T) we will use M = {S^) for the local moment mag- 
netization in this section. We note that M is not the total 
magnetization M = {S~) + {<t~) of the full system of local 
moments and electrons. But due to the self-consistent de- 
termination of {a^){M) ^ n-f — ni described in Sec. [3] M 
is a unique function of M. Thus we can see the local mo- 
ment's magnetization as the order parameter of the whole 
system. 

The free energy is given by its general formula 



Fm{T) ^Um{T)^TSm{T) 



(22) 



=Um{T)+T—Fm{T) . (23) 

Note that the magnetization is considered as a fixed (order- 
) parameter, i.e. an independent thermodynamic variable. 
Using the product rule and subtracting [/a/(0) = ^j\/(0) 
on each side we get 



MAE,Gf}{E)) 



,(M) 



iz^r(-^\E,Gf}{E))- 

iF^J-\E,Gf}{E)) . 



(20) 



t/M(T)-t/M(0) ^ d Fm(T) -Fm jO) 
T2 dT T 

Integration of this equation with the condition 

Fm{T)-Fm{0) 



lim 



T 



'^m(O) 



yields 



where (S*^) '^'^^ is the magnetization of the impurity system. 
The subsystem's Green's function of component ^ can be 
calculated via the projection |fT9l 



Fm{T) =Um{0)-TSm{0)- 
Jo 



(24) 



(25) 



(26) 



y/2 



G(«)(i?) = 



l^irj^^^ -SAE)){GAE)) 



(21) 



' We did not show the calculation of the expectation values 
and the chemical potential for readability. The expectation val- 
ues have to be calculated with the projected functions G^^\E), 
When we choose C = M we get the (local) Green's func- rr'(S) and Fi">(£) while the chemical potential is determined 
tion of the magnetic subsystem G^'^'' (E) and the system of by the full Green's function Go- (E). 
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Since we know the internal energy Um{T) from the pre- 
vious sections we can calculate the free energy if we have 
access to the entropy at T = 0. 

Equation (|26] | allows us to calculate the full temperature 
dependence of the free energy for a special magnetiza- 
tion. When we calculate it for many values of M we get 
an ensemble of functions {Fji/(T)} from which we can 
get the magnetization dependence at a special temperature, 
{Fm{T)} ^ Ft{M) = Uo{M)-T{iT{M) + So{M)). If 
we find its minimum we get the magnetization at the tem- 
perature T: 



dPriM') 



dM' 



(27) 



M 



d{Uo{M') - T(/t(M') + So{M')) 



dM' 



M 



Since this formula is exact, the quality of the results de- 
pend only on the accuracies of the internal energy and the 
entropy . 

4.1 Zero temperature entropy in the concentrated 
lattice We cannot determine the entropy in general so we 
make some approximations to find an expression for Sm (0). 
The zero temperature entropy is defined as 



M 



(28) 



where Fm is the number of ground states according to 
the magnetization M. First we assume that total number 
of states is given by the product of the number of states 
of the local moment and the itinerant electron systerrQ, 
Fm = F^j^jFfj. This is plausible due to the different orders 
of magnitude of the dynamics of the local moments and the 
itinerant electrons. Thus the total entropy is the sum of the 
single entropies 



5M(o) = 5'7(o) + ^t(o) 



(29) 



The assumption of separable electron/spin states is an ansatz 
which would be questionable, e.g., in systems where elec- 
trons are bound to distinct spins. For those cases other, 
more complex, methods to compute the entropy have to 
be found. 

Let us start with the entropy of the local moments. Again 
due to the higher electron dynamics the local moments see 
no single static electron. Thus we can approximate the ac- 
tion of the electrons on the local moments as an effective 
field Bf^ which is equal for all lattice sites. This is a sim- 
ilar situation to the ideal paramagnet except that the field 
depends on electronic properties and the magnetization M 
itself. We do not have to calculate Bfj directly, but we 
can say that each configuration of local moments {mf} 
(m| = —S,...,S) has the same energy as long as it has 



This means that the number of states of one subsystem is in- 
dependent of a special state of the other. 



the same configurational average M' ~ mf . Just as 

well, configurations with different averages do have differ- 
ent energies. At T = only one energy level contributes 
to a thermodynamic average and therefore only states with 
the same configurational AF. So we can conclude that the 
thermodynamic magnetizations equals the configurational 

one, M = A/'. Thus, for a given M, we have to count all 
possible configurations that have the same configurational 
average to get the zero temperature entropy. The identity of 
configurational and thermodynamical average is of course 
not true for T > 0, but in this case we do not need it, be- 
cause we determine the thermodynamical average by the 
minimization of the free energy. For S* = ^ we can give an 
analytical expression of the entropy 



(30) 



Sm{0) n- 
N 



kB ( (i-M)ln(i-M)- 



+ (^ + M)ln(i + M) 



(31) 



and for S' > | we have to count F^^^ numerically. 
The itinerant electrons see more or less a static potential of 
the local moments. That is why we use the expression of 
the Fermi gas 



^^',(0) = - fcs ^ [(1 - (nk.))ln(l - K.))+ 



+ (n-ko-) In (rikcr)] 



(32) 



to evaluate the zero-temperature entropy. Note that the mean 
values (rikcr) are calculated by the full theory of section|3) 
4.2 Zero temperature entropy in the diluted lat- 
tice We need the same simplifications as in the previous 
section. Therefore nothing changes in the formula for the 
electronic entropy (|32] |. except that we get (rikcr) from the 
CPA Green's function. 

For the local moments we define the magnetic subsystem's 
magnetization at T = as M ~ -^jj with N^i = 

xN. When we still assume that the local moments interact 
only with an effective field created by the electrons (cor- 
responding to an ideal paramagnet or independent spins, 
respectively), we do not have to change the general pro- 
cedure to get So{M). Actually the entropy per lattice site 
stays the same as in the concentrated case. Thus we get 



N Nm 
=> 5j\/(0, x) =xSm{0, X = I) 



(33) 



Note that the disorder itself leads to no (magnetization de- 
pendent) contribution to the entropy in this case, but only 
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Figure 1 (color online) The quasi-particle density of states 
of the concentrated system for various magnetizations 
(S^). The QDOS splits into two subband at E ^ -^JS 
and£' w ij(S'+l). For decreasing {S~) spectral weight is 
transferred from upper to the lower subband or vice versa 
(MCDA). The mean-field QDOS (dashed fines) is only tfie 
sfiifted (by z^i J(S'^)) free QDOS. Wfiich magnetization 
(S^) is actually present for a specific parameter set, has to 
be calculated by the minimization of the free energy. Pa- 
rameters: 5 = |, = leV, J = leV, n = 0.05 



the dilution. Thus special disorder effects on the magneti- 
zation, e.g. found in Ref. fSTI . are beyond the scope of this 
paper They could, e.g., play a role in materials which tend 
to form clusters. 

It is a big advantage that we do not need to map the KLM 
onto a Heisenberg model to calculate the local moment's 
properties, especially in the diluted system. This would in- 
volve extra disorder llT0ll22]|23l within the Heisenberg model 
to be calculated separately from the electronic disorder al- 
ready kept within the CPA formalism. Just as the general 
inconsistency of treating one system with two models (KLM, 
Heisenberg) this problem is avoided in our treatment. Since 
our theory is based on translational invariance within each 
component we do not get percolation effects. Those can be 
crucial at low doping if there is a short range interaction be- 
tween the local moments. But even for an interaction range 
slightly larger than next-neighbor, the percolation thresh- 
old is strongly reduced ll24l . Since the interaction between 
local moments in the KLM is assumed to be RKKY-like 
and therefore long ranged, percolation effects should play 
a minor role. 

As mentioned before, in systems which are not transla- 
tional invariant, the entropy had to be calculated more ap- 
propriate, which would be a difficult task indeed. 

5 Results Our main interest is focussed on the depen- 
dence of the Curie temperature on the significant model pa- 
rameters. That's why we first want to discuss these depen- 
dencies in the concentrated {x = 1) system. After that we 
come to the diluted one. Finally we want to use model pa- 



rameters which are suitable for the real system of Gai-^^Mn^jAs 
to compare our theoretical results with those of the experi- 
ment. 

It turns out that in the KLM the term 9a//t(A/) in Eq. (|26] | 
is much smaller than 9a/S't(M). Thus we can neglect this 
in a very good approximation. Since 9a//t(A/) was the 
only term with an explicit temperature dependence we get 
a simplified formula 



T{M) 



dM'UoiM') 



dM'SoiM') 



(34) 



M 



and its inverse function M{T). It means that the magneti- 
zation curve is only dependent on zero temperature quan- 
titiefl Especially we can use ( l34b to calculate the Curie 
temperature, when we set il/ = 0+. 
The lattice is assumed to be simple cubic in this model 
study with a next-neighbor tight-binding dispersion 

W 

e(k) = — (cos(fc2:a) + cos(fcya) + cos(fcza)) (35) 

where W is the free bandwidth. 

5.1 Concentrated System To calculate Tc the quan- 
tity of interest is the internal energy, in our case. Thus we 
should have a look on the quasi-particle density of states 
(QDOS) p„{E) = -l/(7riV)^^ImGka(£;), since it has 
the main influence on U as it can be seen in Eq. (|5]l. In the 
concentrated system the QDOS splits into two subbands at 
different energies (cf. Fig.[T). One is situated at E' « —^S 
(lower subband for J > and upper for J < 0) and the 
other at £' w ■f ("S" + 1). In contrast to the MF approxi- 
mation these subbands remain at this energy range for all 
values of (S*^), only the spectral weight is shifted from the 
upper to the lower subband for decreasing magnetization, 
or vice versa. This means a much lower energy difference 
AU = UiiS"") + d(S'^)) - U{{S'-)) in the MCDA than 
in the MF approximation. Comparing this with Eq. (|34] | 
shows that the MCDA results should give lower Curie tem- 
peratures as the MF approximation. Indeed this is the case, 
e.g. when we plot Tc over J (Fig. |2]i. The MF Curie tem- 
peratures are much higheiQ than those of the MCDA. Ad- 
ditionally we find that ^ while we get saturation 
or only a slight change at large \J\ within the MCDA. 
The next important parameter is the electron density n 
or the hole density p = 2 — 7i in the conduction band. 
In the MFA, due to the shift of the majority spin's QDOS 
by -1^(5'^), the energy difference AU is always positive, 
also resulting in a positive Tc for all n,p (one can also 
see this direct connection in the shape of the Tc -"-curve 



Zero temperature does not mean the absolute ground state 
in our work, since we calculate U for different magnetizations, 
which are not equal in energy 

* Because of the high Tcs of the MFA we usually use a smaller 
I J| for MF results in this work to have a better comparison to the 
MCDA values. 
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Figure 2 (color online) Curie temperature vs. coupling J 
for various electron densities n. In contrast to the mean- 
field result Tc saturates for large \J\ or is growing slowly, 
respectively. At the MF approximation (dashed line, n = 
0.05) Tc increases quadratically with J. Parameters: S ~ 
\,W = leV 
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Figure 4 (color online) Curie temperature vs. free band- 
width W for various couplings J. In the strong coupling 
limit {JS > W) Tc increases approximately linearly with 
W. At lower couplings the effect is inverse and an increase 
of W leads to a reduction of Tc- Parameters: S" = |, 
n = 0.5 
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Figure 3 (color online) Curie temperature vs. electron den- 
sity n for various couplings J. The Curie temperature of 
the MCDA vanishes at half-filling (?i = 1) even for large 
J contrary to the mean-field approximation (dashed line, 
J = 0.05). Parameters: 5 = |, I^ = leV 

in Fig. [3] which resembles the shape of the free QDOS). 
Contrary to that the MCDA has always a vanishing Tc 
at half-filling (cf. Fig. [3]l. Similar results can be found in 
other publications III4ilI0ll25l . This is a very crucial result, 
contrary to a MF approximation which shows a maximum 
Tcs at half-filling (dashed line in Fig. [3]). The reason for 
this behavior is that for n ^ \ the chemical potential /i 
is at the upper band edge of the lower subband. As Fig. 
[T] shows this upper band edge is shifted to lower energies 
for decreasing (5^) within our theory which means a re- 
duction of the internal energy. Thus the energy difference 
dMUo{M) in Eq. (|34] | gets negative resulting in a negative 
Tc, being thus non-physical. 

As the third parameter remains the free bandwidth W. It is 
known 1251 that in the strong coupling limit {JS > W, the 



subbands are split) Tc becomes independent of J (double 
exchange limit). In this case the Curie temperature is pro- 
portional to W as can be seen in Fig. |4] But at lower J 
an inverse effect occurs and an increase of W leads to a 
reduction of the Curie temperature Tc ^ JS/W because 
the bands are not split any more. This is the typical RKKY 
behavior which is valid at low couplings. 
We now can determine the optimal parameters for high 
Curie temperatures in the KLM. Tc increases with the mag- 
nitude of the coupling | J|. But for large | J| saturation of 
Tc occurs and the increase is limited. The electron (hole) 
density should be at quarter-filling for strong couplings or 
lower for a smaller \J\. Finally the free bandwidth should 
be preferably large as long as system stays in the strong 
coupling regime (^ > 1). 

Indeed the findings in this section are similar to studies 
done with a mapping of the KLM onto an Heisenberg model f26l 
I27ll28ll29l . Also in these approaches, which go beyond mean- 
field, a saturation of Tc with increasing \ J\ and a van- 
ishing of ferromagnetism below/at half-filling have been 
found. Thus both methods support each other As stated 
before (Sec. 14.2b the use of the free energy has some ad- 
vantages for computing the magnetization of diluted sys- 
tem. Of course one can get additional information from the 
free energy itself, like the entropy or heat capacity. It also 
can be useful to compare free energies of different mag- 
netic phase to make statements about their existence at fi- 
nite temperatures. This will be left for future work. 

5.2 Diluted system The dependencies of Tc on the 
model parameters in the diluted system stay approximately 
the same as in the concentrated one. But there are some ad- 
ditional effects due to the dilution. 

In diluted systems the carriers are on lattice sites with or 
without magnetic impurities. Due to the interaction with 
the locahzed spins two distinct parts of the QDOS of 
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Figure 5 (color online) The quasi-particle density of states 
of the carriers for various dilutions x. The correlated sub- 
bands around E w ^^JS and E « ^J{S + 1) arise from 
the interaction of the carries with the localized spins Si and 
have a total spectral weight of 2x. Around E « yNM _ q 
a non-correlated sub-band exist resulting from the lattice 
sites without magnetic impurities. Note that the ferromag- 
netically saturated state, which is shown here, is not nec- 
essarily thermodynamic ally stable for all carrier densities 
n and or couplings J. Parameters: S = |, = leV, 
J = leV, Tq^m = 0, n = 0.01, (5") = S 

the carriers arise (Fig. 5). One is located around « (or 
rpNM^ resp.) and results from the non-impurity sites (non- 
correlated part). The two other sub-bands are at w — ^ JS 
and E w ^J{S + 1) at the same positions as the quasi- 
particle sub-bands of the concentrated KLM (correlated 
part)Q. The positions of the sub-bands are not affected by 
the concentration x, but the spectral weights are. In Fig. 
|5]one sees that the correlated subbands have the spectral 
weight 2.T and the non-correlated 2{1 — x). This means, as 
long as the correlated bands are split from the non-correlated 
band, that we get an effective filling n'^^^ = n/xof the cor- 
related band and, e.g., half-filling of the correlated bands 
occurs at a; < n < 2 — x. As seen in the case of the concen- 
trated KLM {x ~ 1) ferromagnetism at half-filling of the 
sub-bands is not thermodynamically stable. Thus there are 
positive Curie temperatures only for n < x (partly filled 
lower correlated subband, cf. Fig. |6) or n > 2 — a; (partly 
filled upper correlated subband) . For x < n < 2 — x 
the chemical potential /i is in the non-correlated band (for 
— at least) andthe lower/upper correlated subband 
is completely filled/empty resulting in a vanishing of fer- 
romagnetism. 

Secondly the maximum Tc at quarter-filling is reduced for 

^ Note that Fig. 5 shows, for graphical reasons, a scenario of 
strong interaction JS/W and all parts of the QDOS belong to 
the itinerant carriers' band (valence band). The correlated sub- 
bands should not be confused with the d-bands of the DMS which 
are taken into account by the local moments in our model. The 
splitting in real DMS would be much smaller. 




n 



Figure 6 (color online) Curie temperature vs. electron den- 
sity n for various dilutions x. A finite positive Curie tem- 
perature only occurs below half-filling of the correlated 
band (n < x). Parameters: S = = leV, J = 0.5eV, 



decreasing x. It can be understood, in the strong coupling 
regime, with the reduction of the effective bandwidth of 
the correlated bandl&J, w ^/xW. Since we know that 
Tc is proportional to W for large \J\ the Curie tempera- 
ture should be Tdn''^) « y/xTdn, x = 1) in the diluted 
system. Note that the strong coupling regime is now de- 
fined by the effective bandwidth, JS > W^^. Thus with 
decreasing x a system can get into the strong coupling limit 
for constant J. 

One of the best known DMS is Gai-^cMn^jAs. To model 
this material we fix some parameters from now on. As 
seen from photoemission experiments |30) the coupling is 
about J sa — leV. Furthermore we approximate the band- 
width from ab-initio calculations OTl to = 4eV. Since 
we use a simple cubic tight-binding dispersion we take 
into account only the main part of the "ab-initio"-DOS be- 
low/around the fermi-level for this estimation. Furthermore 
it is useful to change to the hole picture, where we define 
the hole density p = 2 — n since the bands are completely 
filled for n ~ 2. 

Of course, in a real material other mechanisms besides 
magnetic interactions play a crucial role. For example the 
atomic levels of different atoms will be different in general. 
To take those effects into account we can change the posi- 
tion of the non-correlated band by choosing a specific cen- 
ter of gravity Tq^, which we have introduced in the model 
Hamiltonian (|2]i. Depending on the value of T^'^ the mag- 
netic subbands can be split from the correlated band or can 
lie within it. Especially this determines which kind of holes 
are created by filling the band with carriers. When the up- 
per magnetic subband lies above the non-correlated band, 
holes preferentially occupy magnetic sites for low doping 
leading to a large "magnetic hole density" p'^ (cf. Fig. |7]l. 
This will be reversed when the magnetic subband gets to 
lower energies. The reduction of p^ with increasing 
has two different effects on Tc- For low carrier densities 
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Figure 7 (color online) Electronic QDOS of the diluted 
system for two different centers of gravity of the non- 
correlated band T^^. The total QDOS (black line), the 
projected QDOS of the magnetic sites (red dashed line) 
and the Fermi level (black vertical line) are shown. Top: 

= — 0.5eV. The upper subband stemming from the 
magnetic sites is slightly above the noncorrelated band 
and the lower one is broadened within the noncorrelated 
band. Holes exist above the Fermi level (grey area) al- 
most completely in the magnetic QDOS leading to a ratio 
p'^/x = 0.357 (ratio for all holes p/x = 0.5) . Bottom: 

= 0.5eV. The upper impurities' subband is within the 
noncorrelated band. The number of holes in the magnetic 
QDOS is reduced to p^/x = 0.07 . Parameters: S" = §, 
W = 4eV, J = -leV, a; = 0.05, p = 0.025, {S^) = S, 
T = 0K 



(left of maxima of the parabola-like curves in Fig. |6} a de- 
creasing p^ means a reduction of Tq- On the other hand, 
systems with p > x reach a finite Curie temperature not 
before a certain amount of holes on magnetic sites is trans- 
ferred to nonmagnetic sites (cf. Fig. |8). Figure |9] shows 
clearly the dependence of Tc on the filling of the magnetic 
subband. As long as p^^ /x does not reach a critical value 
a finite Curie temperature exists. For the case that the cor- 
related subband is far above the non-correlated band and a 
large | J| this critical filling would be approximately one. 
If the subband comes near (but not into) the non-correlated 
band hybridization effects with the non-correlated band be- 
come larger and a small amount of spectral weight is trans- 
ferred to the energy region of the non-correlated band. Thus 
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Figure 8 (color online) Curie temperature vs. the center 
of gravity of the noncorrelated band Tq^ for various hole 
densities p = 2 — n. As long as chemical potential is in 
the correlated band Tc is large. When the correlated band 
(and therewith /i) gets into the noncorrelated band (here at 
T^;"^ « (-0.5 . . . 0.5)eV) decreases for low hole den- 
sities. For hole densities p > x finite Curie temperatures 
occur only in the non-split regime. Parameters S — ^, 
W = 4eV, J = -leV, x = 0.05 



the spectral weight of the correlated subband outside the 
non-correlated band decreases below x. Nevertheless the 
complete actual filling in Fig. |9]can be connected to the 
kinks in the p^ — p curves which are at the same posi- 
tion where Tc reaches zero. These kinks disappear when 
the upper correlated subband gets into the non-correlated 
band. This means that a complete filling of the subband 
cannot be achieved even at hole densities p k, x, but at 
much higher hole densities p ^ x (cf. Fig. |7]l. Thus the 
Curie temperature stays above zero for p > x. Even though 
in this case ferromagnetism occurs over a broad range of 
the hole density, the maximum Tc is heavily reduced. A 
very similar picture has been found for dynamical mean- 
field/Monte Carlo calculations lfTSl . 
In fact ab-initio calculations predict that some amount of 
states, arising from the doping with manganese, is slightly 
above the non-correlated band llSTl . This is in agreement 
with measurements of Clio et al. |l8ll9l on the behavior of 
ferromagnetism. They state that for low hole density an in- 
crease of p leads to an increase of Tc- Beyond a specific 
p Tc decreases and finally breaks down as calculated in 
this work in the "split-up" regime. Unfortunately Cho et 
al. cannot give absolute values of p. Why do experiments 
and ab-initio calculations show a maximum at p w a; 
as found by Jungwirth et a/.Q in contrast to our model 
theory? One reason can be that the correlated bands are 
partly within the non-correlated band which allows a finite 
Tc even at p/x « 1 , because the correlated band is not 
completely filled in this case. But more important is that 
the pc?-levels of the impurities are threefold degenerated. 
One hole in these three states mean effectively a "third fiU- 
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Figure 9 (color online) Top: Curie temperatures vs. hole 
density p = 2 — n for different centers of gravity Tq^. 
Bottom: The according hole density on the magnetic sites 
p^/x. Tc is heavily affected by the filling of the magnetic 
subband (cf. Figs. |2l [8]). When the upper correlated band 
is spHt from the non-correlated band (T^^ < -0.25eV), 
p^/x increases rapidly with p. In this case a distinct criti- 
cal ratio can be found where ferromagnetism breaks down 
{Tc ~ 0). If the subband is within the correlated band 
there is ferromagnetism over a broad range of hole densi- 
ties but the maximum Tc is reduced. Parameters: S = ^, 
W = 4eV, J = -leV, 2; = 0.05 



ing" of the correlated band ll32l . Thus even for p ~ x the 
correlated band is not completely filled, but on the left side 
of the maxima of the parabolas in Fig. |9]with an effective 
filling of p'^* w i. This would indeed explain the increase 
of Curie temperatures with increasing hole density and the 
breakdown of ferromagnetism at large p. 
Since we have seen that highest Curie temperatures oc- 
cur in the "split-up" regime and ab-initio calculations also 
predict this scenario we choose a Tq^ where the bands 
are split to investigate the influence of manganese doping 
on Tc- In this case an only partly filled correlated band 
is necessary for finite Tc-s. As the real variation of the 
hole density with doping is not completely clear we pro- 
pose three different types of fillings/compensations in Fig. 
[TOl The first is a constant compensation c resulting in a 
hole density p = x — c. For x < c there is of course no 
ferromagnetism possible since the ordering mechanism is 
carrier mediated. With increasing doping the system has 
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Figure 10 (color online) Curie temperature vs. doping 
rate x. Three curves are plotted for different compensation 
types. A constant hole density (red line, p = 0.02) or com- 
pensation c (black line,p = x — c ~ .T — 0.02) leads to finite 
Tc for X > por X > c, respectively. When p is larger than 
half-filling of the correlated band Tc decreases. A relative 
compensation p = xc (green line, p = .t/2) prefers ferro- 
magnetism at every doping as long as J is large enough (cf. 
Fig.E]) Symbols are from the experiment^. Parameters: 
S=\,W ^ 4eV, J = -leV, T"^^ = -0.5eV 



the optimal p ~ xji where the maximum Tc is reached 
(cf. Fig.|9l). After that further impurities decrease the Curie 
temperature again. The second type of compensation shall 
result in a constant p. Also in this case we need a finite dop- 
ing X to get a positive Tc since x has to be larger than p. 
But with increasing x the Curie temperature changes only 
slightly. Finally there could be a constant relative compen- 
sation p = ex. We have chosen c = i in this work. Thus 
the hole density is for all x near its optimal value and there 
is a large finite Tc- Since there is the same effective fill- 
ing 7 = p/x ~ c and the system is in the strong coupling 
limit, the dominant parameter is the effective bandwidth 
W^^ ~ y/xW. The Curie temperature scales therewith al- 
most only with the effectice bandwidth as it is usual in the 
strong coupling regime (cf. maxima of Tc in Fig.|6]l. This 
results in the curvature Tc ^ \fx- 

Of course these compensation mechanisms are idealized 
and in reality there is probably a combination of these three 
basic types. But it can be seen that the hole compensation 
is indeed a very important effect for the DMS. 

6 Summary We have performed a model calculation 
to describe carrier mediated magnetic ordering of local mo- 
ments. To do that we used a Kondo lattice model with 
ferro- or antiferromagnetic coupling J between the carriers 
and the local moments. The electronic subsystem is treated 
in an equation of motion method to calculate the Green's 
function and the internal energy. From the internal energy 
we get the free energy which can be minimized accord- 
ing to the magnetizion of the local moments. This yields 
the actual magnetization of the system. To include disor- 
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der and dilution of magnetic impurities we have extended 
this method by a coherent potential approximation and a 
proper alloy analogy. 

With this method we investigated the influence of several 
model parameters on the Curie temperature. We found that 
the optimal values are around quarter-filling of the con- 
duction band in the strong coupling \ J\S > W regime. 
Additionally the free bandwidth W should be as large as 
possible as long as the strong coupling limit is fullfilled. 
If W is too large {W > \J\S) Tc decreases. The effect 
of the bandwidth plays also a leading role in diluted sys- 
tems. Since the effective bandwidth of the correlated band, 
which arises from the interaction of carriers and (diluted) 
local moments, is proportional to the square root of the 
doping, « y/xW, increasing x enhances the Curie 
temperature. At half-filling of the correlated band no fer- 
romagnetic ordering exist, which shows the importance of 
hole compensation in diluted magnetic semiconductors. It 
is also necessary for high Tps that the correlated band is 
separated from the non-correlated band and that the chem- 
ical potential is within the correlated band. Finally we have 
compared the calculated values of Tc with experimental 
ones of Gai_j;Mn2: As which are in good agreement to each 
other 

Since we calculate the free energy it is possible to deter- 
mine other quantities like the entropy or heat capacity of a 
system. Also free energies of different (magnetic) phases 
could be compared to see which phase is actually present 
at finite temperatures. We will leave this for future work. 
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